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We have analyzed heat capacity and thermal conductivity measurements of S^RuCU in the nor- 
mal and superconducting state and come to the conclusion that an order parameter with nodal lines 
on the Fermi surface is required to account for the observed low-temperature behavior. A gapped 
order parameter is inconsistent with the reported thermodynamic and transport data. Guided by 
a strongly peaked dynamical susceptibility along the diagonals of the Brillouin zone in neutron- 
scattering data, we suggest a spin-fluctuation mechanism that would favor the pairing state with 
the gap maxima along the zone diagonals (such as for a d xy gap). The most plausible candidates are 
an odd parity, spin-triplet, /-wave pairing state, or an even parity, spin-singlet, d-wave state. Based 
on our analysis of possible pairing functions we propose measurements of the ultrasound attenuation 
and thermal conductivity in the magnetic field to further constrain the list of possible pairing states. 

PACS numbers: 74.25.Fy, 74.25.Bt, 74.25.Ld LA-UR:00-1398 



I. INTRODUCTION 

The search for the superconducting pairing symmetry 
in the layered perovskite material Sr2Ru04 (SrRuO), and 
its attempted theoretical predictions, show remarkable 
parallels to the heavy-fermion superconductor UPt3. In 
both systems, early specific-heat measurements showed a 
large residual value of C /T at low temperatures and were 
interpreted in terms of a superconducting phase with a 
nonunitary p-wave order parameter .EJcl The observation 
of a strong T c suppression with nonmagnetic impurities 
was an additional indication of a superconducting phase 
with an unconventional order parameter.cTu However, 
newer measurements on high-quality single crystals have 
shown that the most likely pairing state in UPt3 is an 
/-wave state, or more precisely a spin-triplet state whose 
orbital basis function belongs to the E 2u representation, 
of the hexagonal crystallographic point group (D6h).E3'EJ 
The experience with UPt3 suggests that the early identi- 
fication of the pairing state, based on low-quality, inho- 
mogeneous samples, is at best inconclusive (for a review 
on UPt3 see, for example, Refs. |l2],[l3]). However, with 
improving sample quality it becomes feasible to identify 
the pairing state by studying transport properties. 

Here we analyze new heat capacity measurements on 
high-quality single crystals of SrRuO, as well as thermal 
conductivity data on dirty samples with a strong Tq-sup-. 
pression, to show that the proposed p-wave model,tirt£l 
A(p^) ~ (jp x + ip y )z, is inconsistent with the available 
data. Our conclusion is that the pairing state in SrRuO, 
most likely, has lines of nodes with gap nodes given by 
the d X y gap function. This can occur in either an /- 
wave state, i.e., a spin-triplet pairing state belonging to 
the E u representation of the tetragonal crystallographic 
point group (D^) or in a d xy singlet state. We argue that 
the /-wave nodal state is consistent with measurements 
of the heat capacity,El thermal conductivity,! 7 ]!^ pene- 
tration depth,£§ Andreev refiection,E3 NMR^Knight 
shift,E3 and /iSR expcriments£3 



Recent band-structure calculations by Mazin and 
Singhcj indicate that there is an increase in the spin sus- 
ceptibility x(q, cj) at four points in the Brillouin zone 
at approximately qo ~ (±27r/3, ±27r/3) that occur due 
to strong nesting effects of quasi-one-dimensional bands 
(£ and £)■ Nesting effects among these bands lead to 
the increased interaction between particles on the Fermi 
surface neac-qoj see Fig. [IJ. In recent neutron-scattering 
experimentsEj the predicted four incommensurate peaks 
near qo were indeed observed thus supporting that nest- 
ing effects near these points are important £3 
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FIG. 1. |£ermi surfaces in the Brillouin zone after Mazin 
and SinghEj The plotted order parameter (proportional to 
opens a gap along (±7T, ±vr) where the incommensurate 
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peaks of the spin susceptibility are observed. The correspond- 
ing quasi-one-dimensional model bands £ at (fc^,, ±2-7r/3) and 
£ at (±2tt/3, k y ) are shown as thin lines. 



In this paper we propose: 1) To identify the regions 
at the Fermi surface near q with the ones that develop 
the largest gap. We use the neutron-scattering data as 
an indication that near the nesting regions the particle- 



1 



particle (or particle-hole) interactions are dominant and 
that these are the regions that would benefit the most 
from opening a superconducting gap. 2) We suggest that 
regardless of the singlet or triplet nature of the pairing 
in SrRuO the gap function should be proportional to a 
d xy harmonics. Such an order parameter would lead to 
lines of nodes along the fc z -axis in the gap and to power- 
law behavior in the thermodynamic and transport prop- 
erties. Line nodes on the Fermi surface lead in clean 
superconductors, and for scattering in the rBjoxii limit, 
to the well-known temperature dependencesc3~E3 of the 
specific heat C ~ T 2 , the nuclear spin relaxation rate 
1/Ti ~ T 3 , the deviation of the penetration depth from 
its zero-temperature value AA ~ T, the thermal con- 
ductivity k ~ T, and the longitudinal sound attenuation 
o:l ~ const., as well as for the transverse attenuation 
a T ~ T 2 . 3) Based on the proposed line nodes in the 
gap we make predictions for ultrasound attenuation and 
thermal conductivity measurements that can further dis- 
tinguish between the remaining possible basis functions. 
We propose complimentary longitudinal and transverse 
attenuation measurements that can help to locate the lo- 
cation of the nodal lines of the order parameter on the 
Fermi surface. Another crucial experiment is the thermal 
conductivity with an in-plane magnetic field. We ex- 
pect the fourfold modulation of the thermal conductivity 
k(6,H) as a function of the angle between the nodes of 
the gap [along the (1,0) and (0,1) direction] and the field 
directions. Thermal conductivity measures the unpaired 
quasiparticle heat transport and is therefore sensitive to 
the angular (field) dependence of the quasiparticle scat- 
tering rate, which "knows" about the angular dependence 
of the gap. We use the analogy with the suggested d- 
wave paring state in high-T c superconductors where this 
fourfold modulation has been observed. l2tEj 



II. MODEL 

The gap function for even parity (spin-singlet) or odd 
parity (spin-triplet) representations is described by an 
order parameter of the form 

A a/3(P/) = A(p/)(iov) a/3 , (singlet) (1) 
Aa/3(P/) = A(p/) • {i(TO- y ) a p , (triplet) (2) 

with o~ a being Pauli matrices. Since nonunitary statesJE@ 
i.e., A x A* ^ 0, have been ruled out by the very 
small residual, value of the specific heat C/T at zero 
temperature]2-3 we restrict our study of spin-triplet states 
to unitary order parameters that factorize into a sin- 
gle spin vector and an orbital amplitude, i.e., A(p/) = 
d A(py), where d is a real unit vector in spin space and 
A(p/) is an odd-parity orbital function. The vector d 
defines the axis along which the Cooper pairs have zero 
spin projection, e.g., if d||z, then Ajf = A^ = and 
A u = A iT = A( P/ ). 



Whether or not spin-orbit coupling is weak or strong 
in Sr2Ru04 has important ramifications for both spin 
and orbital components of the order parameter that are 
allowed by symmetry. While spin-orbit coupling is be- 
lieved to be strong in the heavy-fermion system UPt3 
there are no experimental indications that this is likewise 
true for Sr2Ru04. In the meantime we will use the clas- 
sification of basis functions in terms of irreducible rep- 
resentations of the tetragonal point group (D±h) listed 
in Table |, implying that spin-orbit coupling is strong. 
Since the band-structure palculationsE-3 and de Haas - 
van Alphen measurements^ show very little dispersion 
along k z , we will consider only two-dimensional (2D) ba- 
sis functions on a more or less cylindrical Fermi surface. 
A similar list of possible basis functions was recently com- 
piled by Hasegawa and co-workersE3 for further investiga- 
tions. The listed hybrid state (#3) of the direct product 
B2 g ®E u — E u is a non-trivial realization of the E u repre- 
sentation (also referred to as /-wave state) . So far Knight 
shift data with an in-plane magnetic field H||[100] show 
no change below T c and have been interpreted in terms 
of spin-triplet pairing with the spin vector d locked to 
the crystal c-axis.E3 On the other hand, muon spin rota- 
tion (/LtSR) experiments observed a spontaneous internal 
magnetic field on entering the superconducting state,E2l 
consistent with a time-reversal symmetry breaking state 
belonging to the two-dimensional E u representation. 

TABLE I. 2D polynomial basis functions for the irre- 
ducible representations of D^h of several pairing models (af- 
ter Yip and Garjj23). Notice that B g x E u = E u . The com- 
monly proposed p x + ip y state belongs to the two-dimensional 
E u representation. We present both singlet, d xy , and triplet 
states, B2 g <8> E u , which have lines of nodes, as plausible can- 
didates for Sr2RuC>4. For simplicity we list only the nodal 
angles on the dominating 7 and a Fermi sheets. 
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At this place a caveat is warranted because neither 
Knight shift data at high fields and for a single field orien- 
tation, nor /iSR measurements in impure samples provide 
a clear-cut identification for spin-triplet pairing or bro- 
ken time reversal symmetry states. For example, in UPt3 
early /iSR measurements indicated broken time-reversal 
symmetry in the superconducting phase (probably due 
to impurities), while newer measurements, on very clean 
samples fail to detect any effect at all£3 What makes 
the interpretation of the Knight shift data in SrRuO for 
magnetic fields parallel to the planes even more compli- 
cated, is, that (1) the experiment was not performed in 
the low field limit, but rather deep in the mixed phase, 
H ~ H C 2/2, where contributions from the vortices may 
be important, and (2) nonlocal and surface effects may 
be relevant due to a small Ginzburg-Landau parameter 
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for in-plane currents, 



Ail/en ~ 2.6. 



where 7r,- 



III. THERMODYNAMIC AND TRANSPORT 
PROPERTIES 

We calculate the specific heat and thermal conductiv- 
ity for the order-parameter models listed in Table | and 
fit the results to existing experiments. This way, we can 
determine the model parameters and make predictions 
for sound attenuation measurements. It is important to 
point out that none of the here analyzed transport exper- 
iments can distinguish between a spin-singlet and a spin 
triplet order parameter. Thus we obtain identical results 
for the states #1 and #3. 

The specific heat-pC = TdS/dT, can easily be obtained 
from the entropyOCJ 
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by numerical differentiation. Here /(e) = 1/[1 + 
exp(e/T)] is the Fermi-Dirac function and N(e) = 
— (Nf/ir)lm J dpfg R (pf,e) is the density of states per 
spin with Nf being the normal-state density of states at 
the Fermi surface. 

In the limit of Born (weak) or unitarity (strong) impu- 
rity scattering the in-plane thermal conductpc^ of uni- 
tary spin-triplet superconductors is given byCJ'cil 
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with the unit vector of the Fermi velocity, Vfi, and C R — 
— — -\/|A| 2 — (e R ) 2 . The quasiclassical equilibrium Green 
functions are g R = e R /C R and i R = -A/C R . Within 
the i-matrix approximation for isotropic scattering the 
impurity renormalized quasiparticle energy is e R = e — 
c^pO)- For weak scattering a R mp (e) = (T/tt) J dp f g R , 
and for strong scattering cr R np (e) = — T/(tt J dpjg R ), 
with the normal-state scattering rate T = H/2t. 

In the hydrodynamic regime, cut <C 1, and long wave- 
length limit, ql <C 1, the absorption of ultrasound of 
polarization e pjEqjjagating along direction q is related to 
the viscosity by 



u ,, „ „ „ 
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with the speed of sound c s , the mass density g, and the 
viscosity tensor evaluated at u> — * 0, 
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Here we confine our discussion to order parameters 
with vanishing averages, J dp/ A(p/) = 0, which satisfy 
the gap equation for triplet (singlet) pairing interactions, 



A(p/)= J ^tanh^ |dp^(p / ,p' / )Imf fl (p' / , ( 



(8) 



Note that for spin-singlet pairing all vector functions get 
replaced by the corresponding scalar functions. In the 
weak-coupling spin-fluctuation model the pairing inter- 
action is written as 



V(p f ,p> f )~V*(p f) p' f ) X (pf-p' f ). 

x(q) = xo/[i + C 2 (q-q ) 2 ] . 



(9) 
(10) 
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The detailed form of the effective pairing interaction 
V*(pf,p'f) depends on the form of the spin singlet 
or spin triplet pairing interaction. \q is the static 
spin susceptibility, £ is the antiferromagnetic correla- 
tion length, and the incommensurate wave vectors are 
qo ~ (±27r/3, ±27r/3). The spin-fluctuation scenario pro- 
posed here is similar to the one studied by many, au- 
thors in the conteptt of the heavy-fermion systems,c3 the 
high-T c cupratesxp the quasi-two-dimensional organic 
superconductors,C3 and even SrRuO.E-3£3 In contrast to 
the microscopic model calculations in Refs. [24||45| , we pro- 
pose the existence of either an attractive triplet /-wave 
or singlet d-wave pairing channel in order to describe the 
power-laws observed in thermodynamic and transport co- 
efficients. Our approach is guided by neutron-scattering 
data of the spin susceptibility that can lead to a gap 
function that is gapped along the (±7r, ±7r) directions 
and has nodes along (7r,0) and (0,7r). The immediate 
consequence of the proposed state is that the supercon- 
ducting gap on the holelike (3 band develops nodes at 
(±27r/3,±7r) and (±tt, ±2tt/3). 



IV. RESULTS AND DISCUSSIONS 



In our analysis of the thermodynamic and transport 
properties we make the simplifying assumption that all 
three Fermi surfaces (a, /3, 7) simultaneously go super- 
conducting and can be described by one effective, cylin- 
drical band. At the present time we cannot rule out 
any admixture of the p-wave state #2 to the /-wave 
state #3, since both gap functions belong to the same 
two-dimensional representation E u . However, from a de- 
tailed analysis of the calculated heat capacity we find, 
rather conservatively, that the admixture of a nodeless 
p-wave state has to be lesSpthan 20% to be consistent 
with the experimental C(T).E3 Thus, we neglect the pos- 
sibility of a p-wave admixture to the /-wave gap function 
in the remainder of this work. Impurity calculations for 
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the p-wave state #2 also were performed by Maki and 
PuchkaryoVjCll who reported reasonably good agreement 
between experiment and the calculated p-wave order pa- 
rameter. Very recently, Dahm, Won, and MakO dis- 
carded the nodeless p-wave state and argued in favor of 
/-wave pairing. 

In the temperature range T* <C T <C T c , where 
T* is the characteristic temperature of the impurity 
band width, and in the clean limit, T <C Ao, the eval- 
uation of the entropy and transport coefficients sim- 
plifies significantly. In the presence of line nodes on 
the Fermi surface the density of states is N(e) ~ 
(e/Ao)iV/-. Similarly, we obtain for the Fermi surface av- 
eraged integrand J dp//C(p/,e) ~ — er(e)/Ao, because 
of /C « 7r 4 /(Ree^Ima|^ p )ImC fl , with i R : w e + i0+. 
Thus, the transport coefficients show the usual power 
laws of clean superconductors when using the approxi- 
mate relations for the scattering self-energies, l/2r(e) = 
(-7r) _1 Imcr^ np ~ Te/Ao in the Born limit, or l/2r(e) ~ 
TAo/e for unitarity scattering. 



A. Specific heat 

The states #1 and #3 with line nodes yield C ~ 
NfT 2 1 /aAq, in excellent agreement with experiments, 
while the gapped state #2 disagrees with the data. The 
proposed multiband— |Order-parameter model by Agter- 
berg and co-workersEj, which assumes that only one band 
(7) out of three possible bands goes superconducting at 
T c , fails to describe the low-T dependence (see Fig. ||). 
Our result for the p-wave state #2 is in agreement with 
calculations of the heat capacity by Agterberg.ea In the 
multiband model the density of states (DOS) of the 7 
band is weighted with 57% of the total DOS, while the 
remaining a and [3 bands account for 43% of the to- 
tal DOS. It is the 7 band on which the p-w&ve state 
#2 has been proposed to nucleate. The a and (3 bands 
remain normal. Here fi is the slope parameter of the 
gap function at the nodes, [i = |dA(</>) / Aod(f>\ no d e - In 
our calculations we have used variational basis functions, 
A(p/) — > A(p/)^ r A lg (p/;/i), where the variational func- 
tion Tax belongs to the A\ g representation and remains 
invariant under all group transformations. The slope pa- 
rameter fj, allows us to adjust the opening of the gap 
function at the nodes, which is otherwise not determined 
by symmetry. This enables us to quantitatively describe 
the ground state of the superconducting order parameter 
as probed by low energetic quasiparticles. An approach 
that has been quite successful in describing the low en- 
ergetic quasiparticle excitations in UPt3.EJ 

Assuming that pure SrRuQ has an optimal transition 
temperature of T c q ~ 1.51 K, El we obtain an excellent fit 
for scattering in the Born limit with a scattering phase 
shift <5o — > and a scattering rate T/ttT c q = 0.01. On 
the other hand, resonant scattering (Sq — » 7r/2) with the 
same scattering rate gives a residual value of C/T that is 



too large. If impurity scattering is indeed resonant, then 
a value of T/ttT c0 < 10~ 3 is required to account for the 
lowest measured values of the specific heat. Furthermore, 
it would imply that the optimal transition temperature 
is closer to T c0 ~ 1.48 K. 
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FIG. 2. The specific heat normalized at T c for pairing 
states with line nodes (#1 or #3). A scattering phase 
shift of So — 0° (Born) or So = 90° (resonant), a scatter- 
ing rate F/nT c o — 0.01, and a nodal parameter /1 = 1.5, 
were assumed. For comparison the p-wave state #2 for Born 
(long-dash) and resonant (dot-dash) scattering and the multi- 
band state by Agterberg in the Born limit (cross-dot) are 
shown. The data are from Ref. 

The T c transition of the two components of the triplet 
p-wave order parameter #2, or of the two components 
of the /-wave order parameter #3, is doubly degener- 
ate. Similar to the multicomponent superconducting or- 
der parameter in UPt3 uniaxial strain (pressure) in the 
plane would lifL_the degeneracy of the two-component or- 
der parameter .Ej As a consequence the transition temper- 
ature will split into two. This is a crucial testpefpthe mul- 
ticomponent nature of the order parameterO'EJ Along 
the same line of arguments, a magnetic field in the plane 
should also split T c , as was pointed out in Ref. ^2. 



B. Thermal conductivity 

The in-plane thermal conductivity is isotropic for all 
order parameter models listed in Table ||, assuming a 
cylindrical Fermi surface. In the clean limit, T' « T « 
T c , and neglecting logarithmic corrections, ku ~ T for 
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weak scattering and ~ T 3 for strong scattering. In the 
dirty limit, T <C T* <C T c , the thermal conductivity is 
linear in temperature, Ku ~ T, and independent of the 
scattering strength.. Unfortunately the samples studied 
by Suderow et alx-H exhibit a very strong T c suppression. 
The reported resistive transitions for samples #2 and 
#4, T c e (#2) « 0.81 K and T c e (#4) w 0.58K, occurred 
significantly above the bulk superconducting transitions 
identified by the thermal conductivity, T C K (#2) w 0.60 K 
and T C K (#4) w 0.47 K. Not only does this suggest that 
the samples are in the dirty limit but also that they are 
considerably inhomogeneous. Thus the standard scat- 
tering i-matrix analysis in terms of pointlikc defects in 
the dilute limit will most likely fail to give a quanti- 
tative description. Nevertheless, combining the facts 
of the T c suppression and that the ratios of the resid- 
ual resistivities and the normal-state thermal conductiv- 
ities are related to the scattering rates, r(#41lr(#2) ~ 
M#4)/eo(#2) ~ k n (#2)/k n {#4) « 1-25,EZI we find 
that the normal-state scattering rates are approximately 
given by T(#2)/nT c0 « 0.20 and r(#4)/7rT c0 « 0.25 (see 
Table [il] for the corresponding T c suppression) . 
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FIG. 3. The thermal conductivity normalized at T c 
for the pairing state #1 or #3. The scattering 
phase shifts are <5o = 0°,90°, the scattering rates are 
T/tyTco = 0.01,0.10,0.20,0.25, and fx = 1.5. The data are 
from Ref. [It]. Note that in these dirty samples the resistive 
transition occurs at much higher temperatures (see arrows). 



TABLE II. T c suppression due to the pair-breaking effects 
of nonmagnetic impurities after Abrikosov and Gorkovfij 
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In Figs. H and [| we show the best fits of^^ for sam- 
ples #2 and #4 measured by Suderow et alO. Although 
we cannot obtain a quantitatively good fit for any of the 
pairing models, we are able to ascribe the large residual 
value of k/T to impurity scattering (see Fig. ||) without 
having to invoke a multiband order parameter model (see 
Fig. Q). A surprising result of these fits is that, generally, 
we find better agreement between theory and experiment 
for weak impurity scattering in the Born limit. Very re- 
cently, Tanatar et al.E3 reported measurements of K on 
cleaner crystals (T c w 1.4 K) that arc in good quanti- 
tative agreement with the gapless states Jtl or =#=3 and 
impurity scattering in the unitarity limit O 

For the predicted pairing states #1 or #3, we expect 
to observe a fourfold oscillation of the thermal conduc- 
tivity when a magnetic field is parallel to the layers and 
rotated within the layers. However, the amplitude of the 
oscillations depends on the scattering strength. It is ap- 
preciable for strong scattering (unitarity limit) and very 
small for weak scattering. (Born limit). So far, no oscilla- 
tions have been observed.t§E3 Certainly the experimental 
and theoretical situation remains unresolved and requires 
more study. Indeed such magnetic oscillations have been 
reported in the cuprate YBa2Cu307,ETE3 and are con- 
sidered as additional proof in support of the 4r 2 -j, 2 sym- 
metry of the superconducting state. 
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FIG. 4. The thermal conductivity normalized at T c for the 
multiband model by Agterberg based on the p-wave state 
(#2) and for the same parameters as in Fig. ^| 



C. Sound Attenuation 

The longitudinal (q||e||[100]) and transverse (q| | [100] 
and e\\ [010]) sound attenuations are identical for the pair- 
ing state #2, i.e., a xx (T) / a xx (T c ) = a xy (T)/a xy (T c ). 
This result also was reported in Ref. Whereas for 
pairing states #1 and #3 the longitudinal attenuation 
with q||e||[100] is the same as the transverse attenuation 
rotated by tt/4 with q| | [110] and e||[110]. These rela- 
tions follow directly from Eq. (Q) and are a peculiarity of 
the 2D Fermi surface and the 2D basis functions of the 
order parameters. Inspecting the momentum dependent 
weighting factors in Eq. (Q), 

^L = |cos 2 2,/), 7^ = 1 sin 2 20, (11) 
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it is clear that by rotating the crystal (or the transducer) 
by 7r/4 around the c-axis one simply exchanges these 
functions, w^x n xy ancl j thus swaps the expressions 
for the longitudinal and transverse attenuation. Since 
the integrand /C(p/, e) for the p-wave state (#2) is inde- 
pendent of p/, the longitudinal and transverse attenua- 
tions are identical (within an overall scaling factor due 
to differences in the speed of sound) for arbitrary tem- 
perature and impurity concentration. These predictions 
should be straightforward to check experimentally. In 
Fig. |B| we show the predicted transverse and longitudinal 
sound attenuations for the <i-wave (#1) and /-wave (#3) 
order parameter models. Our results ape similar to the 
ones discussed by Moreno and ColcmanES for the case of 
the d x 2_ y 2-wave gap function in the high-T c cuprates. 
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FIG. 5. The longitudinal and transverse sound attenuation 
normalized at T c for the states #1 or #3, T/nT c o = 0.01, 
So = 0°,90°, and /i = 1.5. In panel (b) the crystal (or de- 
tector) has been rotated by 7r/4 around the c-axis relative to 
the arrangement in panel (a). 



V. CONCLUSIONS 

We have proposed a spin-fluctuation model based on 
the measured spin susceptibility by neutron scattering 
that leads to nodes of the gap function on the Fermi 
surface. We demonstrated that the measured specific 
heat and thermal conductivity are consistent with a spin- 
singlet order parameter (cfcy-wave symmetry belonging 
to B2 g ) or a spin-triplet order parameter (/-wave symme- 
try belonging to E u ), though inconsistent with a gapped 
spin-triplet state (p-wave symmetry belonging to E u ). 
Based on this analysis we proposed sound attenuation 
measurements and thermal conductivity measurements 
in a magnetic field to locate the nodes on the Fermi sur- 
face, as well as measurements of the specific heat sub- 
jected to a uniaxial strain field in the plane in order to 
split the superconducting transition. It is clear that more 
experiments are needed to investigate the nodal regions 



on the Fermi surface and the spin structure of the order 
parameter. 
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